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Abstract 

A small breaking of R-parity reconciles thermal leptogenesis, gravitino dark mat- 
ter and primordial nucleosynthesis. We find that the same breaking relaxes cosmo- 
logical bounds on the axion multiplet. Naturally expected spectra become allowed 
and bounds from late particle decays become so weak that they are superseded by 
bounds from non-thermal axion production. In this sense, the strong CP problem 
serves as an additional motivation for broken R-parity. 



1 Introduction 



A consistent cosmology has to explain the observed matter composition of the Universe. 
Standard thermal leptogenesis [1] provides a simple and elegant explanation for the 
origin of matter. However, thermal leptogenesis cannot be reconciled with an unstable 
gravitino in supergravity [2, 3] unless it is very heavy. Fortunately, the gravitino is 
one of the best-motivated candidates for dark matter. If it is the lightest superparticle, 
the condition that relic gravitinos do not overdose the Universe yields an upper bound 
on the reheating temperature [4]. It is remarkable that a reheating temperature of 
O [lO^^ GeV) can account for the observed matter from leptogenesis and at the same 
time for cold dark matter in the form of thermally produced relic gravitinos with a mass 
of O (100 GeV). Interestingly, the strong CP problem can be solved by the Peccei-Quinn 
mechanism [5, 6] in this cosmological setting. However, the next-to- lightest superparticle 
(NLSP) becomes so long-lived that the strongest conflict between both notions arises 
from NLSP decays that spoil the success of primordial nucleosynthesis [7]. This and 
axion physics lead to tighter constraints on the axion multiplet [8]. 

It has been shown that in the case of small R-parity breaking thermal leptogenesis, 
gravitino dark matter and primordial nucleosynthesis are naturally consistent [9]. Since 
a consistent cosmology should also enable a solution to the strong CP problem, we study 
in this letter how far the restrictions on the axion multiplet are softened, if R-parity is 
broken. 

The small R-parity breaking couplings allow the lightest ordinary superparticle 
(LOSP) to decay into pairs of Standard Model particles. These decays happen instanta- 
neously compared to the Hubble time after LOSP freeze-out or the R-parity-conserving 
decay into the gravitino. Thus they do not endanger the success of Big Bang nucleosyn- 
thesis (BBN). In this way the aforementioned NLSP decay problem is circumvented, 
which at the same time relaxes the constraint on the axino lifetime and thus the allowed 
range of the axino mass and possibly of other parameters. It is allowed to decay right 
before BBN, instead of the requirement with conserved R-parity to decay before NLSP 
freeze-out. Thus its decay temperature can be lowered by three orders of magnitude. In 
addition, saxion decays are allowed to produce NLSPs at any time before BBN. 

The effect on the LOSP decay is the crucial impact of broken R-parity on the con- 
straints on the axion multiplet. All members of the axion multiplet (axion a, axino a 
and saxion (psa^) obtain additional R-parity breaking couplings, but R-parity breaking 
interactions and decays are suppressed by the Peccei-Quinn scale and additionally by the 
small R-parity breaking couplings. Thus they are produced and decay in the same way 
as if R-parity were conserved. One exception to this statement occurs, if the axino is the 
NLSP. In the following we fix the gravitino mass to be of O (100 GeV) and the reheating 
temperature at a rather large value of O (lO^*^ GeV). This is a worst case scenario for 
the axion multiplet. We will comment on the DFSZ [10, 11] and KSVZ [12, 13] invisible 
axion models, especially regarding the possibility of an axino next-to-NLSP. Neverthe- 
less our analysis is general and can be applied to axions in superstring models [14]. Our 
main results are comprised in Table 1. Since, for example, naturally expected spectra 
become allowed, our results serve as an additional motivation for broken R-parity. 
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2 Axino and saxion with conserved R-parity 



Axino Even if axinos do not enter thermal equilibrium after inflation, they are, nev- 
ertheless, regenerated by thermal scatterings and decays in the thermal plasma. The 
resulting density can be estimated in units of today's critical density as [15, 16, 17] 
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= 7.8 X 10- J l^—j— J . (1) 

where and Tr denote the axino mass and the reheating temperature, respectively. 
The axion decay constant fa = fpq/N with /pq denoting the Peccei-Quinn scale if N 
denotes the number of different vacua. For the KSVZ (DFSZ) model = 1 (6). The 
produced LOSP^ and/or gravitino abundances from axino decay are orders of magnitude 
larger than the thermal abundances even if re-annihilation is taken into account [18, 19]. 
The produced LOSP abundance is inconsistent with the scenario, because the NLSP 
decay problem is worsened. Therefore, the decay is demanded to happen before LOSP 
freeze-out, so that one does not have to worry about the produced number of LOSPs 
since they thermalise normally. In this case the lower bound on the axino mass reads 

m,>6xl0 GeV(^l-^j (l^^J [To^^ ) [~^^ ) ' 

where nig denotes the gluino mass, g^, the effective number of relativistic degrees of 
freedom of the Universe and we have used that the freeze-out temperature Ij^^p ~ 
miosp/25. By T-'^'^ we denote the temperature of the Universe after the axino decay. 
Since in supergravity the axino mass is — as the other superparticle masses — generically 
of the order of the gravitino mass [20], i.e., ~ J^T-susy ~ ""^3/2) the mass bound (2) 
can be interpreted as an upper bound on the axion decay constant, i.e., fa < 10^'^ GeV. 
In (2) the decay into a gluino-gluon pair with width [15] 

- 167r3 /2 ml ' ^ ' 

is assumed to be dominant and the strong coupling constant as{fi) = o:s{ma) — 0.1. 
The involved operator exists for any axion model that is able to solve the strong CP 
problem. Investigating other axino decay channels in the KSVZ model, we find that 
they require to be larger than the expected gluino mass, which keeps (3) to be the 
dominant decay channel. Independently of fa the axino must be sufficiently heavier than 
the gluino, which is expected to be among the heavier superparticles due to the running 
of its mass. For instance, with a gluino mass = 1 TeV and the parameter values 
appearing in (2) the axino mass is required to be larger than about 1.35 TeV. 




^The lightest ordinary superparticle (LOSP) is the lightest superparticle of the minimal supersym- 
metric standard model (MSSM). In this sense, gravitino and axino are extraordinary superparticles. 
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The DFSZ superfield is directly coupled to the Higgs superfields. Therefore, the 
DFSZ axino may dominantly decay into a Higgsino-Higgs pair. Using L = j^{H^H^ + 

H^H^) a + h.c, we obtain in the decoupling limit^ 



(4) 



where /x is the usual supersymmetric /i-parameter, tan/5 = Vu/vd is the ratio of the 
Higgs vacuum expectation values and orio and N o^o denote the Higgsino fractions 
of the i-th neutralino. We consider the case of degenerate, maximally mixed Higgsinos, 
which are then decoupled from the gauginos. For sufficiently large tan /? the mixing 
factor in (4) then becomes about 1/2. Thus, the sum of both Higgsinos softens the mass 
bound to 



m^'- > 5 4 X 10^ GcV ( ^" V ( WOGeVy / m,^,, ^ / ,,(T,d 
m~ > 5.4x10 GeY^-^^^j J KWCeV ) 



'dec ^ 

lOO"" 



(5) 



The upper bound on the Peccei-Quinn scale from late axino decay may thus be softened 
to fa ^ 10^^ GeV and for fa < lO^'' GeV we observe that the open decay channel into 
Higgsino-Higgs suffices in the DFSZ model. Other decay channels are subdominant 
compared to (3) and (4). 

Altogether, only spectra with > m-g are allowed in all axion models. This can 
be viewed as a problem, because it requires the axino to be heavier than naturally 
expected. The bound may become softened to nia > mg + ruh in the DFSZ model for 
a small /i-parameter. 

Saxion Due to supersymmetry the couplings of the saxion have the same strength as 
the axino couplings. Thus the saxion is produced as efficiently as the axino. Its dominant 
decay producing MSSM particles is into a pair of gluons with [22] 

2 S 

ygg — '"'sax (f^\ 

^ sax oo_S f2 ^ ' 

or — for the DFSZ model — into a pair of light Higgses with [23, 21] 

W-sax / , 4m^ \ V \ I 



8^ V "^sax/ \faj V^isax/ 

Since in these decays no superparticles are produced, they may happen right before BBN. 
Since the saxion receives its mass from SUSY breaking, one expects msax ~ ""T-susy If 



^The result seems quite insensitive to the possibility of the Guidice-Masiero mechanism to generate 
the /i-term, but this might be different for the saxion decay into two light Higgses [21]. For the neglected 
phase space factor we find {(1 + — Xh)^{l + + Xh)^{{i — x^f — xl)}^^^ with Xx = m^o/m~ and 
Xh = mh/m~. 
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the saxion is heavy enough to produce superparticle pairs, its decay could lead to the 
same worsening of the NLSP decay problem as the axino decay. In this situation the 
lower bound on the saxion mass from the decay into a pair of gluons becomes in the end 



mosp \i / fa Y f9*{Ts^^)^~' 



If the DFSZ saxion decays dominantly into Higgs pairs, we can derive an upper bound 
on the saxion mass 

4m2 \ ^ /iQii GeVy / ^ \4/l02GeV\V 100 ~' 



TTls 



< 740 GeV 1 



IJ \ fa J UoOGevJ V "inisp J W{Ti^ 



(9) 

We see that for fa ^ 10^*^ GeV the DFSZ saxion mass is only required to allow for the 
decay into a pair of light Higgses. 

In addition, the saxion may decay as well into two axions with [24] 

2 S 

-paa ^ 'sax (^r.\ 

647r /2 ' 

where the self-coupling x can be of order 1. The produced axions represent a form of 
dark radiation, i.e., decoupled, relativistic particles not present in the Standard Model. 
During BBN the energy density of dark radiation pdr is constrained to be less than the 
energy density of one additional neutrino species [25], which translates into 



Pdr 



<0.14. (11) 

BBN 



Here psM denotes the energy density as expected from the Standard Model and 

PdrlBBN - ^aa PsaxlTd- \ J J ' ^ ' 

where g^:s counts the relativistic degrees of freedom determining the entropy density s 
and g*s{T^^^) = g*{T^^^). Written as bound on the branching ratio Baa of the saxion 
into two axions (11) reads'^ 

Baa<0A{l + 50.xy[^-^j[^^^) fcJl^^J • ^''^ 

In this inequality we have approximated the decay width of the saxion as Fgax — Tfax + 
Fg^^. Thus our conclusion should hold qualitatively for any axion model. Then the 
branching ratio reduces to 

■paa 2 
R ZIM - Z (14) 



Our expression differs from the literature in the numerical value and the dependence on . This is 
because we have taken into account properly the scaling of pdr oc gti^T'^. The simplifying assumption 
of entropy conservation (g.s ~ const.), which is usually made, does indeed not hold during these stages 
of the Universe. 



5 



The value for fa appearing in (13) corresponds to the upper bound on fa from axino 
decay and the axion energy density ila (see below). Since in our scenario the reheating 
temperature is fixed at rather large values, all members of the axion multiplet enter 
thermal equilibrium after inflation for such small values of fa- Therefore the saxion 
yield l^ax cannot be smaller than the equilibrium value Kjll ~ 1.21 x 10~^. The ap- 
pearing values for fa and rrisax are chosen to show the worst situation in the considered 
scenario. Like a smaller fa or a larger ?nsaxi also a larger x leads to an earlier decay, 
which corresponds to a smaller Osax at its decay. Thus there is a self-curing effect for 
large x. The bound (13) represents indeed an implicit equation for the self-coupling x. 
Evaluating it for x it turns out that there is no constraint on x at all in the scenario 
under consideration. 

We point out that the absence of any bound on x is due to the expectation mgax ~ 
msusy and the restriction of fa to small values appropriate for the considered scenario. 
Furthermore, Ygax could be much larger than Kjax, if the saxion is produced from coherent 
oscillations after inflation as discussed below. 

The saxion might decay before LOSP freeze-out, if the self-coupling is strong enough. 
Neglecting conservatively the saxion decay into gluons and Higgses we find 

^ V "^sax J VlQio GeVy' I102 GevH 61.75 J ' ^ ' 

In this situation there is no additional constraint on the saxion mass. Note that the 
decay into axions does not produce a significant amount of entropy. 

If a new best-fit value demands additional radiation energy in the Universe [26, 27, 28, 
29] , which is often parameterised by a change of the effective neutrino degrees of freedom 
AA^eff > 0, we can determine parameter values from (13), such that the additional 
energy is formed by axions from saxion decay. For fa < lO^*' GeV and m-sax ^ 10^ GeV 
the maximal AN^s is 0.6. However, when these requirements are relaxed also larger 
AA'efj are possible. For example, we obtain AA'gfj ~ 1 with a rather small saxion mass 
nT-sax = 10 GeV, fa = 10^° GeV and a self-coupling x = 0.1. 



3 R-parity violating case 

Axino Since in the case of broken R-parity the NLSP decay problem is absent, the ax- 
ino may decay right before BBN. Produced superparticles decay promptly into particles 
of the Standard Model, which thermalise normally. Then the lower bound (2) becomes 

,„,>89Gev(l-^) (l5T4v)'(iSv)'' 

where we still assume the decay channel a ^ gg to dominate and take T^^^ as lower 
bound on the temperature after the particle decay. The bound T^^^ = 4 MeV, corre- 
sponding to a lifetime Tmax ~ 0.05 s, considers the neutrino energy density during the 
neutrino thermalisation process [30, 31, 32], if the particle dominated the energy density 
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of the Universe before its decay. This is hkely to happen in the considered scenario (see 
below). Weaker bounds T^^f^ ~ 0.7 MeV arise from BBN calculations [33, 34]. As in 
the following we omit the dependence on g^{T~'^'^ = = 10.75. We see that due to 

-^min ^ -^n?sp axino mass bound becomes so weak that masses much smaller than the 
expected gluino mass would become allowed. 

In this situation the axino decays always depend on which channels are kinematically 
open and thus, in principle, on the full spectrum. This is a qualitative difference to the 
R-parity conserving case. Particularly interesting is the possibility of an axino next-to- 
NLSP, which we will assume in the following. Then its dominant decay is fixed into the 
NLSP, in other words, into the LOSP. From (16) we see that the NLSP were allowed 
to be a gluino. A recent analysis reports m-g > 322 GeV as experimental lower limit for 
such a long-lived gluino [35]. From (5) we see that the DFSZ model would in addition 
allow for a Higgsino NLSP, if the corresponding decay channel were kinematically open. 

Since the lightest neutralino is likely one of the lightest superparticles in the spec- 
trum, one interesting decay is into neutralino and photon with [15] 

where C„^o^ = {CaBBl cos 0m/)A^^oro > while N is the bino fraction of the i-th neu- 
tralino and Qw denotes the weak mixing angle. We take the electromagnetic coupling 
constant aem(^) = Oem("ia) ^ 1/128. The axion to two B bosons coupling CaBB 
varies for different implementations of different axion models.^ For simplicity we set 
CaBB = N^ogo = 1. We find 

,„,>45Gev(l-^) (iSt4v)'(4w)' 

and conclude that a neutralino with a substantial bino component can be the NLSP, if 
the axino is the next-to-NLSP. Thus a superparticle spectrum with nig > > m^o is 
now easily possible without any BBN conflict. 

If the axino decay is fixed into a sneutrino, it decays via an intermediate neutralino 
into a photon and a sneutrino-neutrino pair, i.e., a — )• 'jx^* ~^ ^Vv. This process is 
suppressed compared to (17) by an additional power of ttcm and further factors depending 
on the neutralino composition and the exact spectrum. Thus in case of a sneutrino NLSP 
an axino next-to-NLSP might be possible only for parameter values at the boundaries 
of the allowed region and a tuned spectrum. Consequently, this situation is disfavoured. 

For other sfermion NLSPs the situation depends much more strongly on the axion 
model. In the DFSZ model we estimate the kinematically unsuppressed decay width of 

*For example, in the DFSZ model with (d'^, e) unification, CaBB = 8/3. In the KSVZ model, for dif- 
ferent electromagnetic charges of the heavy quark eg = 0, —1/3, 2/3, CaBB = 0, 2/3, 8/3, respectively. 
Below the QCD scale, CaBB is reduced by 1.92 [36]. 



1287r3 fl 



m 



2 
a 



(17) 
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the axino into a stop-top pair using the Lagrangian of [37] as 



where rrit denotes the fermion mass and = l/(tan^/3 + 1). This decay arises from 
a dimension-four operator that becomes important at low decay temperatures and ac- 
cordingly small masses. If only this channel were open, the resulting lower mass bound 
on the axino would become 

The lower bound would practically be given by fa and the requirement to have the 
channel kinematically unsuppressed, since a long-lived stop with mass below 249 GeV 
is excluded according to the CDF experiment [38] . It might well happen that the decay 
of the next-to-NLSP (a) into the NLSP {t ) and its superpartner (t) is kinematically 
forbidden. The axino would decay violating R-parity. This is the same if the a is the 
NLSP. We comment on this excluded case below. 

The decay width into other sfermions is given by (19), if mt is replaced by the 
corresponding fermion mass and — ?• X^ = tan^/3/(l -|- tan^ /3) if appropriate. For 
example, the decay width into stau and tau is (19) with the replacements nit — t- m-r ~ 
1.78 GeV and X^ — ?• X^. The corresponding lower bound on the axino mass becomes 

Thus, for fa below 10^^ GeV also this channel is restricted by the requirement of kine- 
matic accessibility only. Since nir/m:^ <^ 1, this is not a strong constraint, if the stau is 
not much heavier than the gravitino. Consequently, if the stau is the NLSP, the DFSZ 
axino may be the next-to-NLSP. 

Since the other leptons are lighter, in these cases the bound becomes tighter. While 
the smuon stays possible, the selectron would require > 1 TeV even for parameter 
values at the boundaries. These considerations expand to the quarks lighter than the top. 
The superpartners of the bottom, charm, and strange quarks are possibilities. For the 
down squark the Peccei-Quinn scale needs to be close to the lower limit. The situation 
for the up squark is even worse than the one for the selectron. 

In the KSVZ model axino-sfermion-fermion interactions are one-loop-suppressed in 
the low-energy effective theory [39, 40], which weakens axino to sfermion- fermion decays 
substantially relative to the DFSZ case. The tree-level decay via intermediate neutralino 
into other sfermion-fermion pairs is comparable to the decay into sneutrino-neutrino. 
Consequently, the KSVZ axino is disfavoured, if the axino decay is fixed into a sfermion- 
fermion pair. 

If the axino were the NLSP, i.e., for a spectrum with rn-iosp > fria > it would 

decay either into the gravitino or via R-parity violation. Both decay modes are strongly 



8 



suppressed by the Planck scale or by the Peccei-Quinn scale and the R-parity violating 
coupling, respectively. Therefore, the lifetime of the axino would always become much 
larger than the time of BBN. This is independent of whether R-parity is broken by 
bilinear or trilinear couplings [41, 42, 43]. Since ila ^ 1; such a late decay would spoil 
the predictions of BBN. We conclude that an axino NLSP stays excluded in the R-parity 
violating scenario. 



Saxion If R-parity is violated also superparticle pairs from saxion decays are harmless. 
If its dominant decay channel is into a pair of (massless) gluons, the most severe bound 
on the saxion mass is found to be similar to (16). So for /„ ~ 10^" GeV the saxion mass 
is practically not constrained in the R-parity violating case. However, considerations 
concerning the self-coupling x are not affected by R-parity violation. 

If late-decaying particles enter thermal equilibrium after inflation, they likely dom- 
inate the energy density of the Universe at their decay, since the energy density of 
matter pmat grows relative to the energy density of radiation prad as the scale factor a 
of the Universe, i.e., Pmat/Prad oc a. If a particle dominates the energy density of the 
Universe, it might produce considerable entropy at its decay, which could dilute cos- 
mological abundances to too small values. Axino and saxion enter thermal equilibrium 
at Tr ~ lO^*' GeV if fa < 10^^'^ GeV. The dilution factor due to their decay can be 
estimated to be [44, 8] 

fa \ flTeV\-^ fg^T'''^)^-^ 

^a/sax 



assuming the decay into a gluino-gluon pair (3) and two gluons (6), respectively, to be 
dominant. The numerical difference between axino and saxion is small. Here, the upper 
bound on the dilution factor A > 1 arises from leptogenesis itself and is A™^^ ~ 10^. 
Thus the occurring dilution is so small that there arises no constraint from entropy 
production by the thermally produced axion multiplet. 

The saxion may also be produced by coherent oscillations around its potential mini- 
mum after inflation. The temperature at the onset of saxion oscillations is 

^ 2.2 X IQio GeV ( " ( ^^^1 ' . (23) 
VlTeV; WiT°i^)J 

If the reheating temperature Tr < T°^^ as in the standard scenario, the produced saxion 
abundance is [24] 

osc 

r sax — Tp^ 



Mpi 



4.2 X 10-9 GeV ( ^—^y^) ( -rj^ ) ( ^ ) , (24) 



2 X 109 GeV J VlO^o GeV J \ fa 
where cpl.^^ denotes the initial amplitude of the oscillations and where we have assumed 

2^sax'Asax- 



the simplest saxion potential, V = hm^ 
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The dilution factor due to saxion decays A ~ 0.75 Tg^^^/^sa^ with T^^-^ = |^ 



,osc 



IS 



given by 

Since the dilution factor is by definition > 1, values for A smaller than one lead to the 
standard scenario with A = 1. In this discussion we fix the reheating temperature at its 
lower boundary from thermal leptogenesis, = 2 x 10^ GeV [45] or 2A x 10^ GeV 

with A > 1. The tightest bound on the initial amplitude in the standard scenario is 
found for the maximal axion decay constant fa = 10^^ GeV and a small saxion mass, 
for concreteness mgax = 10^ GeV. It is (pl^^ ^ 7 x 10^'^ GeV. In comparison, the loosest 
bound is found for the minimal axion decay constant /a = 6 x 10^ GeV and a rather 
large saxion mass, for concreteness mgax = 1 TeV. It is (pl^^ < 1.4 x 10^^ GeV. Thus 
both bounds are far below the Planck scale. Bounds on the initial amplitude of saxion 
oscillations are not affected by R-parity violation. Conservatively allowing for smaller 
saxion masses as well, we summarise these bounds in Table 1 as (p].^^ < (lO^^-lO"*^^) GeV. 

In the case of considerable entropy production, for a certain A the reheating tem- 
perature has to become larger than the temperature at the onset of saxion oscillations, 
i.e., Tr > T°i;i. Then Tr in (24) has to be replaced by T°^^ from (23). In this case 
larger initial amplitudes may be allowed, while at the same time larger axion decay 
constants become allowed. The tightest bound on the initial amplitude in the scenario 
with A = A™^^ = 10^ is found for the maximal axion decay constant /a = 4 x 10^^ GeV 
and the minimal saxion mass msax ~ 300 GeV from the decay into a gluon pair. It 
is (?l)gax ^ 4 X 10^^ GeV. In comparison, the loosest bound is found for the minimal 
axion decay constant /a = 6 x 10^ GeV and a rather large saxion mass, for concreteness 
"^sax = 1 TeV. It is (pl.^^^ ^ 5 x 10^^ GeV. We summarise these bounds, that are still far 
below the Planck scale, in Table 1 as a range (pl^^ ^ 5 x (lO^^'-lO^^) GeV. 



4 Constraints from the axion 

First, the lower limit on the axion decay constant from axion physics [46], /a ^ 6 x 
10^ GeV, is not changed by the additional R-parity violating interactions. Furthermore, 
the axion still decays harmlessly into two photons with a lifetime much larger than 
the age of the Universe. Due to its tiny mass, the thermal relic density of axions is 
negligible. In contrast, a large axion density may be produced non-thermally by vacuum 
misalignment and topological defects. 

The density produced by vacuum misalignment [47, 48, 49] is usually given as 

where oq comprises model-dependent factors. In the presence of gravitino dark matter 
we require Qa/^DM = r <^ 1, which gives an upper bound 

fa < 10^° GeV for ao = N = 1 , r = 0.04 . (27) 
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Note that this bound becomes tighter for larger N. 

Topological defects occur, if the Peccei-Quinn symmetry is restored after inflation, 
i.e., Tr > TpQ ~ /pQ = Nfa- This leads to the formation of disastrous domain walls [50], 
if > 1, and cosmic strings [51]. The abundance of axions from cosmic strings Q^^^ 
exceeds that from vacuum misalignment [52, 53]. Since here the reheating temperature 
is fixed at about 2A x 10^ GeV, the situation is particularly interesting. There are two 
possible cases: i) topological defects are not created after inflation, so we do not have 
to care about them or ii) they occur and we have to take them into consideration. To 
avoid topological defects completely, 

/a>2xl09GeV^, (28) 

which could lead, depending on N, to a stronger lower bound on fa than given above. 
In the case without entropy production (A = 1), this favours models with > 4 such as 
the DFSZ model. If (28) is violated, N = 1— fulfilled by the KSVZ model— stih avoids 
domain walls and the axion density from strings, ~ 10 x ^i™^*^ [52, 53], gives a tighter 
upper bound than (27), 

fa < 1.3 X 10^ GeV for Oq = iV = A = 1 , r = 0.04 , (29) 

where Cq comprises model-dependent factors at the production of axions from cosmic 
strings. Combining the above considerations, only the small interval (1.3-2) x 10^ GeV 
for fa might be excluded. In this sense, the allowed band for fa is not changed. 

With entropy production (A > 1) the reheating temperature is raised to compensate 
the dilution of the baryon asymmetry, so the Peccei-Quinn symmetry becomes restored 
for a larger range of values of fa- We consider entropy production after the QCD phase 
transition by late particle decay and estimate the axion abundance from cosmic strings 
as n^"^' ~ 10 X i7™'VA with Q™^^ as in (26). Then the upper bound (29) on is softened 
by a factor of A^/'', because the axions are diluted, while the amount produced remains 
the same. The maximal A™*^^ ~ 10^ corresponds to /a < 4 x 10^^ GeV, so in this 
situation the Peccei-Quinn symmetry is probably restored, because Tr ~ 2 x 10^'^ GeV. 
We therefore list this upper bound on fa in Table 1. Larger values are possible, if the 
symmetry is not restored and if one accepts the axino to be heavier than the gluino. 

Also the upper bound on fa from vacuum misalignment is affected by A > 1 if the 
entropy is produced after the QCD phase transition. In this case the Universe is likely 
to be dominated by matter at the onset of axion oscillations at T"^'^ ~ 1 GeV and the 
upper bound on fa becomes [54] 

( fa y < (^\ /4MeV 
V 1014 GeV y ~ V0.02/ V rdcc 

where T'^'^'^ denotes as before the temperature of the Universe after the entropy-producing 
particle decayed. In this case the axion decay constant is constrained more strongly by 
O^*"" and too late axino and saxion decays. Altogether, a large dilution factor opens up 
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— at least for = 1 — more parameter space, but the situation depends on the time of 
entropy production. 

We would like to point out the appealing possibility that the saxion decay gives 
rise to the maximal dilution after the QCD phase transition. For the harmless value 
fa = 10^^ GeV this is the case for an initial amplitude of saxion oscillations equal to the 
geometric mean of the two scales involved, i.e., (pl.^^-^ ~ VfaMp, and a quite small saxion 
mass msax ~ 10 GeV [8]. However, it is understood that an adjusted initial amplitude 
might produce the desired amount of entropy for any of the allowed parameter values. 

The given bounds from the different production mechanisms of axions refer in each 
case to "standard values" in parameter space. They can be relaxed or circumvented in 
"non-standard" scenarios. Constructing models that realise a small is beyond the 
scope of this letter. 

5 Conclusions and outlook 

We have presented how the constraints on the axion multiplet in scenarios of thermal 
leptogenesis with gravitino dark matter are relaxed, if R-parity is broken. Our results are 
comprised in Table 1. Obviously, the most important findings also hold in less restricted 
scenarios. 

The saxion mass becomes practically unconstrained. An axino NLSP stays excluded, 
but the axino may be anywhere else in the superparticle mass spectrum as long as a decay 
channel into an ordinary superparticle of the MSSM is kinematically open and allowed 
at tree-level in the low-energy effective theory. This does not hold if the only possible 
decay is into a light fermion (up-quark, electron, neutrino) and its superpartner. In 
turn, for the DFSZ axion model, usually considered superparticle spectra are possible, 
including spectra with stau NLSP. The sufhcient condition is that the axino decay into 
the NLSP is not too strongly kinematically suppressed. In the KSVZ model, decays into 
fermion-sfermion pairs are loop-suppressed. However, in both models an open axino 
decay channel into a neutralino suffices. As a consequence, it will be interesting to 
construct concrete models with the naturally expected axino and saxion mass of order 
"T-susyj possibly along the lines of existing models like [21, 55, 56]. 

Constraints on the axion decay constant from axino and saxion decay are softened. 
Interestingly, the DFSZ model is — depending on the /i-parameter — less restricted from 
particle decay already in the standard scenario. For example, the saxion mass could 
be unconstrained, if the Higgs boson is lighter than the NLSP. Since in the considered 
setting the reheating temperature is fixed at a large value, the constraints from axion 
overproduction on the decay constant are particularly interesting and depend on the 
axion model. However, they are not affected by R-parity violation. Therefore we have 
shown how late-time entropy production might soften the upper bounds on the axion 
decay constant and the initial amplitude of saxion oscillations after inflation. 

There is no constraint on the saxion-axion-axion self-coupling as long as fa < 
W^^ GeV and msax > 100 GeV. A large self-coupling of order one can remove con- 
straints on the saxion mass, especially in the standard scenario. We remark that for 
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[GeV] 


standard 






fa 


<10io 


<10io 


< 4 X 10^2 


^a 


> ^9 

(or > + TUh in DFSZ) 


> "i"nlsp" 


> max[m"nisp", 

30o(ioi/gcv)'] 


^sax 


>760(T0i4iv)' 

(or > 2mh in DFSZ) 

2 

^ [5 (iQiu^Gcv) '2mnlsp] 


- ^ (lO« Gev) 


- ^ (lOi'> Gev) 


'/'sax 




< 10l3_iQl5 


< 5 X (IQi^-lQiS) 



Table 1: Parameter constraints for the different scenarios (standard: R-parity conserved 
(A = 1), R-parity violated (A = 1) and violated R-parity with the maximal entropy 
dilution at the right time). Units are GeV where not written explicitly. Here, we assume 
the self-coupling x ^ 1. Only if a mass bound is sensitive to the actual value of fa, its 
dependence is given. By "nlsp" we indicate that the bound actually does not hold for 
any possible NLSP and depends on the axion model. The upper bound on /„ in the 
standard scenario is boldface to indicate that it arises from and the (KSVZ) axino 
decay. In the ^ case it stems from il^^^ only. If, furthermore, matter dominates at T°^^ 
it stems from ^l^^- 

suitable values of the self-coupling and other model parameters one can obtain any de- 
sired amount of additional radiation energy in the Universe formed by axions. It will be 
interesting to identify (more) models realising a particular self-coupling and/or a small 
axion density Qa ^ ^^dm- 

Furthermore, it might be interesting to investigate if other scenarios of dark and 
visible matter enable a solution of the strong CP problem. We conclude that broken 
R-parity does not only solve the NLSP decay problem but also makes it easier to solve 
the strong CP problem. Our results are an additional motivation for scenarios of thermal 
leptogenesis with gravitino dark matter and broken R-parity. 
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